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ABSTRACT 


This report describes the techniques used to determine the 
trajectories and orientations of three sounding rockets instrumented 
to study the aurora. The radar plot board data were fitted to a 
near earth expansion of the force of gravity to determine the 
trajectory. Only onboard magnetometer data were used to determine 
the attitude of the payload with respect to the earth's magnetic 
field. Computer programs in the FORTRAN language are available 
which generate the trajectory and attitude parameters. 
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Chapter I 
Introduction 


A 

We have launched three Nike-Tomahawk sounding rockets from 
Fort Churchill, Canada, to study aurorae. In order to determine pitch 
angle information about energetic electrons and for analysis of data 
from other experiments on board it is necessary to be able to describe 
the position and attitude of the payload. The payloads are launched 
with a Nike booster which falls away when spent and a Tomahawk. second 
stage which burns out before atmospheric exit and remains attached to 
the payload. Flight duration and apogee are typically 500 seconds and 
250 km respectively. To achieve stability the vehicles are spun at about 
7 rps during burning. After burn-out they are despun to about 1 rps to 
facilitate collection of angular information in the data, Radar plots 
furnished by the Churchill Research Range were used to determine vehicle 
position. On board magnetometers which measured the component of the 
earth's magnetic field parallel to their orientation were used to determine 
payload attitude. 


Chapter II 
Vehicle Position 

The effects of coriolis and centripetal accelerations upon sounding 

rockets launched from Churchill show up most dramatically in reducing the 

eastward distance traveled by several kilometers (the earth rotates under 

the payload) # The North-South location of the impact point is virtually 

unaffected. In practice one is most concerned with the effect upon the 
*NASA 18.63 UE, 18.64 UE and 18.65 UE. 
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altitude where the effect of centripetal acceleration is less than 0,1%. 
However, since most vehicles are launched eastward the coriolis accelera- 
tion is in the opposite sense to the centripetal acceleration and can be of 
nearly the same magnitude so one can safely ignore these effects. Because 
the exact impact point is generally not needed and because radar data 
frequently is not good enough to justify further precision we will adopt 
a coordinate system assuming a flat, nonrotating earth, with; positive z 
representing attitude and positive x representing eastward direction. 
Normally a flat earth assumption would have g , the acceleration of 
gravity, independent of z . However over the range of z for the sound- 
ing rocket this represents an appreciable error so we will use an expansion 
of the potential energy, V , for the Inverse square gravitations,! force, 


( 2 - 1 ) 


V = - 


GM m 
e 


- GM m 

e 

R +z 
e 


where G is the gravitational constant, M £ is the mass of the ear tH, 

m is the mass of the payload, r is radial distance from center of the 

earth and R is the radius of the earth at Churchill, 
e 


Expanding ( 2 - 1 ) in z/R g gives 


V = 


- GM m 
e 

R 


-1 


( 1 +z/R ) 


- GM m 
e 


[l-z/R e + (z/R £ ) 2 - (z/R £ ) 3 + ...] 


GM m 
e 

R 


GM m GM m _ 

. e e 2 

+ z 3 z 

R R 


neglecting the cubic and higher terms. 
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(2-3) 


The payload kinetic energy can be written, 

• 2 • 2 *2 

T = 1/2 m (x + y + z ) . 


(2.4) 


Hence the Lagrangian is , 

• 2 *2 *2 GM m 

L - T - V = 1/2 m (x + y + z ) + — — 

R 


GM m GM m 

e , e 

o 2 "f" o 


R 


R 


Since x and y do not appear in L the corresponding 
velocities are constant above the atmosphere. The differential equation 
in z is 


(2-5) 


• 1 

z 




+ 2 g. 



where 


is the acceleration of gravity on surface pf earth, 


(2-6) 


GM 

e 



The general solution to (2-5) is 


(2-7) 


z = B 1 e at + B 2 e _0tt + B 3 


where a and B^ are given by, 

(2-8) a = (2g o /R e ) 1/2 B 3 = R e /2 . 

Using (2-6) and the value of g measured at Churchill of 
2 

g = 981.761 cm/sec one can solve for the value of R to use. These 

e 

values are in Table I. 


J 
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Table I 

R„ = 6377.0 km 

Ej 

= 3188.5 km 

a 2 = 0.307907 x 10 -5 sec" 2 

a = 0.175473 x 10 -2 sec" 1 

The problem remaining is to use the measured z^ and t 
from the radar plot to least squares fit for the initial conditions 
B^ and B 2 . This analysis was quite satisfactory for 18:63 and 
18:64, but for 18:65 we obtained a better fit to the points using 
slightly different values for Rg and fitting for also. 

Apogee time, t , and height, z^ , can be found by differ 
entiating (2-7) 

(2-9) fc A - J CB 2 /B 1 ) 1/2 

(2-10) z a~ - 2 (B 1 B 2 ) 1/2 + B 3 . 

The values for liftoff time, t. , z. and the four coefficients 

A A 

in (2-7) for t and t^ measured from liftoff are given in Table 11, 

Table II 


Vehicle 

Liftoff 

fc A 

z A [km] 

B^[km] 

B 2 [km] 

B 3 [km] 

atsec" 1 ] 

18:63 

21 Mar ’68 
0601:32.7 

241.09 

241.76 

-965.129 

-2249.25 

3188 o 5 

0.17547x10 

18:64 

14 Jan '70 
0405:30.0 

250.24 

259.74 

-943.961 

-2271.7 

3188.5 

0.17547X10 

18:65 

17 Jan '70 
0303:20.0 

248.16 

252.48 

-779.05 

-1989.9 

2742.7 

0.18894x10 
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One can use (2-7) to determine v by differentiation. Figures 

z 

1-3 are plots of the altitude and v versus Universal Time in minutes. 

z 

The x and y components of velocity are given in Table ITT. 


Table III 


Vehicle 

v (East) [km/sec] 

X 

v y (North) [km/sec] 

v xy [km/sec] 

18:63 

0.051 

- 0.101 

0.113 

18:64 

0.203 

- 0.005 

0.206 

18:65 

0.276 

0.168 

0.322 


It is interesting to compare the altitude at which the payloads 
were inverted by the atmospheric drag on the fins and the altitude at 
which the electronics began experiencing sustained high voltage breakdowns 
in Table IV. 



Table IV 


Vehicle 

Turnover altitude [km] 

Breakdown altitude 

18:63 

72 

83 

18:64 

65 

81 

18:65 

73 

81 


Because the radar did not track the payloads throughout the flights 
there may be several kilometer uncertainties in z below about 120 km on the 
downward leg of the flight. The uncertainty on the upward leg is on the order 

of 100 meters. It is possible to obtain an exact analytical solution to this 
problem for t as a function of z , but it would be very difficult to use 
least squares technique to obtain the initial conditions. 
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Chapter III 
Payload Attitude 

After despin the payload is sufficiently above the atmosphere to 

enable one to neglect torques due to atmospheric friction. In the absence 

of net torques L , the total angular momentum vector, is constant in an 

“> f 

inertial frame. We will then choose the direction of L to be the + z 
direction of the space axes for our description of payload attitude. This 
space set of axes is not the system used to describe position of the payload. 
More aspect information than supplied by the magnetometers is necessary to 
relate the two systems. 

Barring the unfortunate and rare case where t and B would be 
parallel or antiparallel we will use the direction of B to provide the 

other direction necessary for the space axes. Assuming that over the altitude 

and temporal range of interest B is nearly constant in direction we define 

the x direction by specifying that B lie in the x' - z' plane and that 

B x , , the component of B parallel to x T , be negative. In a simplified 

case where L is in the local vertical direction and the magnetic declination 

is zero this coordinate system would have the x' axis pointing south 

(equatorward from Churchill) and the y* axis pointing eastward because the 

magnetic field is in the northward direction. Figure 4 represents this 

simplified case. In this coordinate system we can describe the magnetic field 

as 

(3-1) B = B (t) b = B ,(t) i' + B ,(t) k' = B (t) [sin 3 i' + cos 3 k'] 

OX z o 

where the magnitude of B , B q (t) , does have the altitude dependence, 

^ A ^ 

through t , and 0 is the angle between B and k* , the unit vector along L 
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By definition of i* cos 3 < 0 and for a typical L direction sin 3 < 0 
so 3 would lie in the third quadrant. 

The body coordinate system, x" y" z" , is used to describe 
position and orientation within the payload. We choose the z" axis 
to be a principal axis and assume that it coincides with the symmetry axis 
whcih describes the axial dimension of the payload parallel to the geometric 
center line and in the direction of the nose of the payload. The x" and 
y" axes are also principal axes and are assumed along the directions 
specified by the magnetometers as in Figure 5. The origins of the two 
systems coincide at the center of gravity of the payload. 

For time-independent moments of inertia (the slowly despinning 
case will be considered in an appendix) and for a rigid body rotating about 
its center of gravity with body-fixed axes coinciding with the principal 
axes Eulers equations are 


(3-2) 


1 ii ti W ii = M tiW ii (I it it - I it it 

x x x y y y"y" z z 


I II || GO tt GO ||G0 II (1 || || — 1 II || 
y y y z x z z z x 


1 it n GO n — (0 ||(0 || (I || || — I || || 
Z Z Z X y XX y"y" 


) 

) 

) • 


Because of the symmetry about the z" axis we can define transverse, 
I,j, , and axial, 1^ , moments of inertia 


I 

I 


ii 


ii 



I 


T 


(3-3) 
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For our payloads 1^ » I (1^ ^ 200 1^) . Equations (3-2) now can 


be written 


(3-4) 


h v 

X T “y" 

V 


“y"V (I I - V 

- W (I T - I A ) 


= 0 


From the last of equations (3-4) we have 


(3-5) 

w z „ “ 0 ; 

I - I 

to „ = constant . 
z 

Define 

T A 

ft = 0 ) „ ( — ) , then equations (3-4) become 

2 X T 

• 

(i) = (i) ,.ft 

x y 

(3-6) 

• 

CO fi “* CO 

y x 

1 A V = ° * 



The first two equations can 

be combined to give 


• 

00 f, — (a) ,,ft — 
x y 

0 

(3-7) 

• 



0) „ + 03 „ft = 

y x 

0 . 


Multiplying the first equation by eo „ and the second by U) „ and 

x y 


summing them gives 


("US'* 
\ - w / 


fil ..(l) .. + 03 ..0) .. = 0 

x" x" y" y" 


Equation (3-8) can be integrated to give 

(3-9) 

therefore from (3-5) and (3-9) we have, 


2 . 2 „ _ 2 
“x" + “y" = constant = W T 
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(3-10) 




+ 00 



constant, 


the angular velocity vector has constant magnitude. 

Assuming constant fi one can differentiate the first of equations 
(3-6) and substitute in the second one to give 

^2 


(3-11) 


co n — oo nfi — — co i,fi 
x" y x 


Equation (3-11) describes simple harmonic motion for co^,, . For typical 

initial conditions let co „ = 0 and co „ = co using (3-9) . Then we can 

x y i 

write for co ,, and co „ 
x y 

to „ = co sin fit 

(3-12) 

COy, = C 0 T cos fit . 

We can also write co and L in the body frame 


(3-13) 


co = co „ i" + co „ j " + w „ k" 

x y z 

A A A 

= co T (sin fit i" + cos fit j") + co z ,,k" 

L = T-u 

A A A 

= I co (sin fit i"+ cos fit j") + I co ,, k" 

1 1 A Z 


From Figure 6 we can obtain the angle 
z" axis, and £ , the angle between co and the 


(3-14) 


tan 0 = 


I T ai T 

W 


co„ 


co z „ 


0 between L and the 
z" axis. 


(3-15) 


tan £ 
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Since 1^ > I , W lies between L and the z" axis. The angles 
9 and £ are constant because the quantities on the right hand side of 
equations (3-14) and (3-15) are constant. 

(3-16) 0 = e = 0 0 = 8 q , £ = £ q . 


We use the Euler angles to refer the body system to the space 
fixed system. Note that the angle 9 is also the Euler angle 0 (see 
Figure 7). The transformation between the two sets of axes is given by^ 

• • 

0 ) ,, = <j) sin 0 sin ip + 8 cos ip 

(3-17) w „ = | sin 0 cos ip - 0 sin ip 

• • 

« = <j> cos 8 + ip 

Using (3-16) these equations become 

= <p sin 8 sin ^ 

(3-18) i0y M = | sin 0 cos ip 

• • 

w „ = cj) cos 9 + ip 

From (3-9) and (3-18) we have 

2 2 2 *2 2 
(3-19) u „ + w „ = oo = <p sin 0 . 

x y i 


Using (3-16) and (3-19) we can define the precession frequency, 


/.\ 


P 


9 “J 


(3-20) 


<b = constant = to 

P 


From the last of equations (3-18) and (3-5) and (3-20) we can 

define the spin frequency, w , 

s 
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(3-21) ip = constant = to . 

s 

Differentiating (3-18) and substituting them into the first 
equation of (3-4) gives 


(3-22) I T # sin 0 cos ip = (4> sin 0 cos i p) (<£ cos 0 + l p) ( I ^ - 1^) 


which can be solved for <£ in terms of ^ 
(3-23) 


• V 

^ (Irp-Ij COS 0 


t i ) 


T A 


COS 0 I 


and we have the ratio of spin frequency to precession frequency 


(3-24) 


01 / 00 = (— - 1) cos 0 . 

s p 'I 


Using (3-16), (3-20) and (3-21) we can write the Euler angles 
in terms of their time derivatives 

0 = 0 

o 

(3-25) <j> = oi t + <p 

p o 

ip = 0) t + ip 
r s r o 

We are now prepared to describe the way the magnetic field, B , 
which is fixed in the space frame will be seen by the magnetometers on 
board the payload - the body frame. 

The magnetometers measure the component of magnetic field 
parallel to their orientation. They output a bias voltage, CB , of about 
2,5 V plus a voltage linearly dependent upon the magnitude of the component 
of B parallel to their orientation. The linear coefficient, CA , is 
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about 4.0V/Gauss. Let n\| be a unit vector in body frame in the direction 
measured by the i magnetometer, then the output voltage, , has the form 

(3-26) V ± = CA i B • n” + CB i . 

■> A 

Therefore if B is antiparallel to n^ the output is less than CB^ and 
vice-versa. The analysis used here will assume the form of (3-26) for V 
and neglect non-linear terms. 

From (3-26), if B and n'^ were measured in the same frame, the 
output voltage would be a constant term modulated by a cosine term, but alas 
all is not so simple! We need to know the components of B in the body frame 
to be able to compute the dot product required in (3-26). This transformation, 
using (3-1) , is given by 



(3-28) A = 


cos 4> cos <f> - cos 9 sin cp sin ip 

cos ip sin tf> + cos 8 cos (p sin ip 

- sin ip cos (j) - cos 8 sin 4) cos ip 

- sin sin (ft + cos 0 cos j) cos \p 


\ sin 9 sin <p - sin 8 cos 4) 


cos 9 



Equation (3-27) then becomes 


= B „ i" + B „ j" + B „ k" 
x y J z 


(3-29) 


A*B 
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where 


B^,, = B Q (t) [sin 3 (cos ip cos (p - cos 0 sin <p sin ip) + cos 3 

sin ij> sin 0] 

(3-30) = B Q (t) [sin 3 ( - sin ip cos <p - cos 0 sin (p cos i p) + cos 3 

cos ip sin 0] 

B = B (t) [ s * n 3 sin 0 sin <f> + cos 3 cos 0] . 
z" o 


Now one can compute the dot products for the three possible 
magnetometers measuring the x" , y" , z" components of the field. For 

A A 

the x" component, n.„ in (3-26), is i" and 


(3-31) B" • n „ = B (t) [sin 3 cos ip cos (p - sin 3 cos 0 sin <j> sin ip + 

X o 

+ cos 3 sin ip sin 0] 


For the y" component, n^,, is j" 


and 


(3-32) B" • n^,, = B o (t) [- sin 3 sin ip cos (p - sin 3 cos 0 sin (j) cos \p + 

+ cos 3 cos ip sin 0] 


For the z" component n.„ is k" and 


(3-33) B" • n z „ = [sin 3 sin 0 sin (p + cos 3 cos 0] . 


To simplify these equations we define 


(3-34) 


a = 
b - 
c = 
d = 
e = 
f = 


sin 3 sin 0 
cos 3 cos 0 
sin 3 cos 0 
cos 3 sin 0 
d - c sin <p 
sin 3 cos <p 
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We shall restrict ourselves to the case where L is upward rather 
than downward over Chur chill. This is the condition before despin and it is 
improbable that despinning could invert the payload angular momentum vector. 

For 18:64 and 18:65 the z" axis magnetometers indicate that L remained nearly 
vertical. This assumption puts 3 in the third quadrant and makes the first 
four quantities defined in (3-34) negative. 

The resulting equations are 

->■ A 

(3-35) B" • n x „ = B Q (t) [e sin \p + f cos 

(3-36) B" • n^„ = B Q (t) [e cos \{j - f sin \Jj] 

A 

(3-37) B" • n z „ = B Q (t) [a sin <p + b] . 

We choose t in equation (3-25) to be zero at a time when the 
payload (and z" axis) is most antiparallel to B (when z” is nearest -B) . 
This means that equation (3-37) is at a minimum value indicating (a < 0) that 
sin <J> is unity, hence 


(j)(t = 0) = 4> 0 = tt/2 

(3-38) e Q = e(t = 0) = d - c 

f = f(t = 0) - 0 . 

Because (p varies much slower than ip one can also require that 
t = 0 be chosen at a time where V^„ , the voltage from y" ~ axis magnetometer 
( the RAM-) is less than CB^.. - the bias value. This means that (3-36) is at a 
minimum value at t = 0 . Recalling that f = 0 determines to be 0 or 

ir depending upon the sign of e Q . Because c and d are typically both 
negative we can define a positive quantity, g , 


(3-39) 


g = c/d = tan 8/tan 0 . 
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If g > 1 , then e >0 and [Li = n and If g < 1 , e <0 

and tJj q = 0 . Since 6 is in the third quadrant (3-39) implies that if 

8 - tt > 0 , vp Q = tt whereas if g-TT<0,ip o = O. When 8 - it > 0 the 

payload precession cone does not include B and when 8 “ tt < 0 the 

->■ 

precession cone does include B . Hence, if the payload precesses around 
B , i(Jq = 0 and if the payload does not precess around B , ip Q = tt . Figures 
8 and 9 show the two cases. 

We will now outline the procedure for determining the two constant 
angles 0 and 8 . We have measured the times and voltages at local maxima 
and minima on the y" - axis magnetometer and x" - axis magnetometer for 
18:64 and 18:65 respectively. 18:63 had only the y" - axis magnetometer. 

For 18:64 and 18:65 we made detailed, simultaneous measurements from all three 
magnetometers during the several second time interval where the payloads turned 
over and all three magnetometers sampled a full range of values. Taking ad- 
vantage of the fact that the magnitude of the field was nearly constant over 
this short time interval we were able to make self-consistent checks on the 
calibrated values of the CA's and CB's . 

We also picked several times encompassing the complete altitude 
range where we simultaneously measured the voltage output from all three 
magnetometers. Using the calibrations we converted these voltages into values 
for the magnetic field components parallel to the three axes. We then summed 
the components to obtain the field magnitudes. 

(3-40) B o (t) = [B x „ 2 (t) + B^, 2 (t) + B z „ 2 (t)] 1/2 , 


over a range of altitudes. We found that we could model the altitude dependence 
as 


(3-41) 


B Q (t) = B Q (z(t)) 


B 

oo 


(Rg+ztt )) 3 


as would be expected for a dipole field model. 
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The next step is to look at the components of B which are 
parallel and perpendicular to the payload spin axis, z" . Let a be 
the instantaneous angle between the z" axis and B , then we define 


(3-42) 


B^ = | B | sin a 

B 1 1 e | B | cos a 
| B | E (B* + B^ ) 1/2 = B Q (t) . 


Since we measure x" or y" magnetometer voltages only at their 
local maxima or minima the corresponding y" or x" magnetometers are 
perpendicular to Bj_ and measure no field at these times. Hence, by measur- 
ing the voltages at local maxima and minima we know that they correspond to 
the instantaneous value of B^ . The import of this is that in general the 
dot product in (3-26) for the x" magnetometer is 


n i, 
x 


= B „ 
x 


which is only some part of 

we have B • n „ 

x 


Bj- , but at the special times we measured 

= B „ = B. = B (t) sin a . 
x J_ o 



To generalize let m refer to the indices x" or y" when y" 
or x" respectively are perpendicular to the field, then (3-26) becomes 


(3-43) 


V = CA B (t) sin a + CB . 
m mo m 


Equation (3-43) has only sin a unknown. It can be solved for 
sin a and then for cos a , 


(3-44) 


cos a = ± 


1/2 


1 - 


<y c v 

(CA m B o (t#J 
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We will choose the negative sign in (3-44) when we can tell from 

the envelope of the values of V that the z" axis is still above the 

m 

magnetic horizon. This will be discussed in further detail later. 

Of course, cos a is very easy to compute from the output of the 
z" - axis magnetometer, 

(3-45) V „ = CA „ B (t) cos a + CB ,, 

z z o z 

V - C V 

cos a = p-yv 

CA Z " B 0 (t> 

here there is no ambiguity about the sign of cos a . 

Comparison of (3-26) , (3-37) and (3-45) shows that 


(3-46) 


cos a = b + a sin (j> . 


By inspection of the envelope curve for a node or by finding 
minima in B^ one can determine the Universal Time , t Q , where the z" 
axis is closest to - B and the voltage is a local minimum. The time to 
the next similar node gives , the precession period. One can then 
write <f> 


(3-47) 


<P 



t + d> 
Y o 


* * 

where t = t - t and d> = n/2 and t is Universal Time, 
o o 

Armed with an expression for 4> one can perform a least squares 
fit to equation (3-46) to determine the constants a and b . Using the 
first two equations of (3-34) we found two pairs of 3 and 6 (a < 0 , b < 0) 
One pair, 3^ and 8 ^ , corresponds to a precession cone not including B 
and the other pair, 82 an d ^ oes correspond to z" - axis precession 


around B . 
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Vehicle 18:63 had only a y" - axis magnetometer which was shifted 
upwards in CBy, and consequently for much of the flight local maxima of 
V„ were greater than the telemetry voltage limit. We also suffered some 
telemetry dropouts which unfortunately coincided in time with the local 
maxima and minima of Vy, further hampering data reduction. To fill the 
gaps where Vy, could not be measured at local minima or maxima we 
decided to use the values of local maxima or minima of V ,, which were re- 

y 

constructed from the slope of V ,, evaluated at CBy, . We had about 90 
directly measured values and about 50 indirectly measured values of Vy, 
with some overlapping to check the accuracy of the slope reconstruction 
method. 

Using (3-26) and (3-36) we have 


Vy, = CAy, B Q (t) [e cos ^ - f sin ip] + CB „ . 


(3-48) 

The expression in brackets can be written 


(3-49) 


[e cos ip - f sin ^] = h sin (ip Q - ip) 


2 2 2 

where h = e + f and tan ip Q = e/f . Taking advantage of the fact that 
B Q (t) , e and f are nearly constant over a few spin periods one can compute 
the derivativeof (3-48) , 


3V „ 


(3-50) 


9tT = CA y" B 0 ( t ) h COS W 0 " W 


= - CAy, B Q (t) h \p cos (i|i Q - ip) . 


But we measured the slope where V „ = CB „ which meant that ip - ip = mr 

y y ' T 0 

and consequently cos (ip Q - ip) = ±1 . This means that we can solve (3-50) 
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for h in terms of known quantities (ip is the spin frequency which for 
18:63 could be directly measured to a few percent), 


(3-51) 




/ CA „ B o (t) . 


But at a time ± T g / 4 from where V^,, = CB^„ the angle ip will change to 
where sin (ip Q - ip) = ± 1 and we can substitute the value of h from 
(3-51) into (3-49) - (3-48) to get 


(3-52) 



/ ip + 


CB „ 

y 


CB „ 

y 


Therefore by measuring the slope and spin rate we were able to 
reconstruct values to give the correct envelope voltages to use in determin- 
ing 0 and 0 . This method was checked in several overlap regions by 
measuring both the slope and local maxima and minima and the results agreed 
to within 4 %. 

Normally the envelope of local maxima-minima voltages will exhibit 
nodes with the frequency of the precession of z" about L . These nodes 
correspond to times where the z" - axis is nearest to - B . For 18:63 
secondary nodes betweeen precession period nodes indicated that the z" - axis 
had dropped below the magnetic horizon giving the envelope primary nodes 
when (p = tt/2 + 2mr and secondary nodes when <J> = 3tt/2 + 2niT (see figure 10). 
This was very fortunate giving a built-in calibration of the product 
CA^ B Q (t) in equation (3-44) at the times when the payload was perpendicular 
to B . Equation (3-46) had /a/>/b / and we varied the amount of time spent 
below the magnetic horizon to give the best fit to the measured envelope. 
Figure 11 indicates that the fit was very good except when a was near 90® 
when the telemetry problem was most severe. 
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By measuring times of adjacent magnetometer maxima and minima 
one can determine the spin frequency, to g to only about ± 4% because (3-48) 
for example is not a pure sine wave and the time between adjacent maxima 
reflects variations in parameters other than if) alone. To improve this 
measurement we counted the number of oscillations, n Q , in an approximate 
precession period (to the nearest complete cycle). We then accurately 
(± 5 ms) measured the time elapse, At , during the n Q oscillations. The 
hypothesis was then made that to g had one of the three values given by 


(3-53) 



27T(n o -l) 

At 

2fTn 

o 

At 

2TT(n +1) 
o 

At 


A computer program was written to make the final determination of 

the correct set of 0,6 and to g . The procedure involved a double loop 

which tested the six possible combinations of the to g 's and (6 , 0) 's in 

equation (3-35) or (3-36) to predict the times and voltages of the maxima 

and minima. For each of the three possible to 's both (6 1 , 6.) and 

©2 ♦ ^ 2 ^ were tested. The procedure was very sensitive in that the five 

incorrect possibilities had obvious phase shifts from the observed maxima- 

minima during precession whereas the one correct combination predicted the 

maxima-minima times to within the measurement error. The results had the 

pair (3-, , 6.) with to = to and the pair (3 9 , 0 9 ) had to = to 
XX s s z z ss 

*r O 

In order to specify the attitude of the payload with respect to B 
one needs to determine 0 , the coning half-angle, 3 , the angle between 
L and B , <p , the Euler precession angle and ip , the Euler spin angle. 
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4) is found from the primary nodes in the envelope curves. The time between 
them gives , the precession period. (j> is found using (3-47). 0 and 3 

are found by fitting equations (3-46) for a and b and then simultaneously 
solving the first two equations of (3-34) for pairs of (9 , 3). The choice 
of which pair of. (0,3) and which u) g to use is then determined by testing 
the possible cases to minimize the phase and amplitude errors between results 
from equations (3-35) or (3-36) and the measured values. 

Tables V and VI give the frequency-period results and the angular 
* results respectively. 


Table V 


Vehicle 

T(j) (first node) 

TP (Precession period) 

TS (spin period) 

18:63 

0604:34.38 

131.5 sec 

1.93976sec 

18:64 

0409:38.516 

185.751 

1.05540 

18:65 

0306:26.245 

199.104 

1.06209 



Table VI 




average angle 

e , 3 

Vehicle 

-e- 

o 

O 

CD 

3 between z" and -B 

uncertainty 

18:63 

tt/2 0 70.2° 

205.6° 72° 

± 6° 

18:64 

tt/2 tt , 20.50 

203.4 31° 

± 1° 

18:65 

tt/2 tt 9.023 

°193,6° 16° 

± 2° 
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APPENDIX I 

Evaluation .of Euler Angles for Despinning Symmetric Payload 

Vehicles 18:64 UE and 18:65 UE carried experiments which required 
the deployment of long booms. After the yo-yo despin the motors which 
deployed these booms were switched on. Deployment time was about 100 sec. 
The booms extended radially, in opposite directions, from the spin axis. 
Their contribution to the transverse moment of inertia, 1^, was negligible, 
and the vehicle was assumed to remain axially symmetric. The effect 
upon the axial moment of inertia, I , was significant and resulted 
in further decreasing the vehicle-spin rate. To make complete use of 
data from the boom experiments and other experiments during the deploy- 
ment it was necessary to determine the Euler angles for this time period. 

The total time derivative of angular momentum, L , with respect 

< 1 >■ 

to inertial coordinates for a varying moment of inertia, I , and a 

3 

torque, N , is given by 

-+ • • 

(A-l) N = — = I -00 + [ I ] • cu + 0 ) x £ 

at 

where for a symmetrical payload changing only the axial moment of inertia. 


I 



n n 

w W 

A A 

j" 3 i x o 

A A 

k" k" l. J 
A 


0 


0 
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0 


0 0 

0 0 

0 k"k"I A J 


We assume the external torque, N = 0 . The component equations 
analogous to (3-4) are 


I 03 + 03 ,,1 ,03 1,- 03 t ,I_a) „ = 0 

T x y A z z T y ' 


(A- 2) 


I l“y" + “z" W' 


“AV ’ 0 


I a“z" 


Vz" 


+ 03 


,,Im^3 ,, — 0) ,,X_03 |, — 0 
i T y y T x 


From the last of (A- 2) we have 


(A-3) ^z" + V°z" = 0 ’ Vz" = constant = a • 

The magnitude of the total angular momentum, |L| , is a constant. 
Because the component of L along z" is constant the Euler angle 0 
is defined by 


(A- 4) 


sinG 


Vz" 



constant 


Hence, for our case, 0 , the coning angle is constant during the 


long boom deploy despin 
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The first two equations of (A-2) can be manipulated as in (3-7) to 
(3-9) to give. 


(A-5) 


2 , 2 2 
ui „ + a) ,, = constant = oj m 

x y ' T 


However because w „ changes the angular velocity vector doesn’t 
have constant magnitude and the angle £ increases as w M decreases. 


(A-6) 


tan 8 = 


I t oj t I T w T 

*aV = ~ 


UL 


tan e = 


V 


Using equations (3-17), (3-19) and(A-5) we have <f> = constants , 


and 


(A- 7) 


u) „(t) = iJj(t)+<j)cos9 . 
z 


We define 


(A- 8) 


w = i cos0 = constant 
o T 


then using (A-7) , (A-8) and (A-3) 


(A- 9) ijj(t) = w z „(t) - oj q = — - oj q . 


<}> , (}) and 0 can be determined from the rigid body portion of the 


flight. 
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By choosing an appropriate model for I (t) we can integrate (A-9), 

A 

and compute the remaining Euler angle ip . 

The axial moment of inertia can be made from a constant component, 

A , plus a contribution due to the deploying booms, Ig(t) • If the booms 
deploy at a constant speed ^Boo^*-) will have quadratic time dependence. 


(A- 10) 


V £ > ■ B(t - C D y 


I A (t) - A+B(t-t D )' 


where t Q is the time the booms begin deploying. We know that the value 

2 

of A is 2.196 kg-m . Using (A-10) we can integrate (A-9) from t^ to 

V 


(A- 11) ^(t 2 )-ip(t 1 ) = a 


dt 


[A+B(t-t D )*] ** a, o (t 2“ t l ) 


^(tJ-iKO = 


/"AB L 


tan 


-1 ‘VV'® -1 ‘V t D >> ®' 

A tan A J - WV 


We drop the suffix on t 1 and let t 2 = t f , the time the long booms 
complete deployment. 


^ ( t) = —2— tan 
v^AB 


(t-t D )/AB _ 

T - 0) t + Ip 

A J o Y o 


(A- 12) 


where ip = ip(t c ) - ip. + w t- 
o f T f of 


and iK = tan 

f /AB 


-i rvv^J 
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Actually the argument of the inverse tangent is related to the boom 
moment of inertia 


(A-13) 


ARG 


4 


(t-t D )- B 


i 


V') 

A 


We can determine ip(t^) by using (3-25) which is appropriate for 
the part of the flight where I is constant. 

(A-14) iKt f ) = u> s (t f -t o ) + it 


Two more constants, a and B , remain to be evaluated to determine 
(A-12). From magnetometer data we can measure the observed spin period, 
T„(t) . It depends upon cp and 6 as well as ip because during most of 

b 

the flight ip is constant whereas deviates as the vehicle approaches 

its nearest position to the magnetic field. However for times > 0.05 t 

P 

from one of these magnetometer envelope nodes we can write 


(A- 15) 


T s (t W | 

<P 


then using (A-9) we have 


(A-16) 


T_(t) 

5 ' ' 


Ztt 


a-u I. (t) 
o A 


Using (A-10) and making an expansion for a»U) A 

o 


(A-17) 


T .(t) = . 

b a -Ui A 


2TTA 


(a-w A) 
o 


r ....a) a i 

L 1 + 7^Ta)\ 


9 


one obtains 
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Equation (A-17) is now in a form which can be least squares fitted 
by 

(A-18) T g (t) = C 1 +C 2 (t-t D ) 2 . 

Using the values of and determined by the fitting 

we have 


(A-19) 


2uA . 
a = — — + a) A 
o 


(A-20) 


B = 


2tt 


C 0 (a-a) A) 
£ o 

L(a-a)A) + 


The constant B is proportional to the rate at which the booms 
deploy. We found that the spin rate indicated that they deployed at 
two rates, a slow initial rate and then at a higher rate until they were 
completely deployed. Consequently equation (A-12) had to be used twice 
with separate values of a and B . 
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